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Abstract 

The elliptic flow for Z <2 particles in heavy ion collisions at energies from several tens to several 

> . 

tJ" , hundreds MeV per nucleon is investigated by means of transport model, i-e. a new version of the 

O 

I Improved Quantum Molecular Dynamics model (ImQMD05). In this model, a complete Skyrme 

■ potential energy density functional is employed. The influence of different effective interactions and 

medium corrections of nucleon-nucleon cross sections on the elliptic flow are studied. Our results 
show that a soft nuclear equation of state and incident energy dependent in-medium nucleon- 
^ ' nucleon cross sections are required for describing the excitation function of the elliptic flow at 



intermediate energies. The size dependence of transition energies for the elliptic flow at intermediate 



^ \ energies is also studied. The system size dependence of transition energies fits a power of system 

d \ size with a exponent of 0.223. 
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One of the main goal for the research area of heavy ion conisions(HICs) at in- 
termediate energies is to extract more accurate information on the nuclear equation 
of state(EoS). Considerable progress has been made recently in determining the equa- 
tion of state of nuclear matter from heavy-ion reaction data A promi- 
nent role among available observables is played by collective flow. A lot of theoreti- 
cal and experimental efforts on the study of the collective flow in HICs have been paid 
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, Q 18, Q, 3, 21 , Q, Q The elliptic flow has proven 
to be one of the more fruitful probes for extracting the EoS and the dynamics of heavy ion 
collisions. The parameter of elliptic flow is quantified by the second order Fourier coefficient 

2 2 

V2 = (cos 20) = ( p2+p2 ) from the azimuthal distribution of detected particles at mid-rapidity 
as follows 

dN , 

— = Po(l + 2t;icos0 + 2t;2cos2</)) (1) 

(Jjip 

where (f) is the azimuthal angle of the emitted particle momentum relative to the x axis. 
Positive values for (cos 20) reflect a preferential in-plane emission, and negative values for 
(cos 20) reflect a preferential out-of-plane emission. The change-of-sign recently observed at 



ultra-relativistic energies has received particu 



ar interest as it reflects the increasing pressure 



buildup in the non-isotropic collision zone[2J]. Recently, the excitation function of elliptic 
flow parameters at energies from Fermi energy to relativistic energy regime for ^'^'^ Au+^'^'^ Au 
has been measured by FOPI, INDRA, ALADIN Collabrations0,Q| and the transition energy 
from positive to negative elliptic flow was conflrmed, which is around 100 MeV per nucleon. 
The elliptical flow parameters (cos 20) at energies from tens to hundreds MeV per nucleon 
are determined by the complex interplay among expansion, rotation and the shadowing of 
spectators. Both the mean fleld and two-body collision parts play important role at this 
energy region. The mean fleld plays dominant role at low energies and then gradually the 
two body collision part becomes dominant with energy increasing. Thus, a detailed study 
on the excitation function of elliptical flow at this energy region can provide more useful 
information on the nucleon-nucleon interaction related to the equation of state of nuclear 
matter and the medium correction of nucleon-nucleon cross sections. The transition energy of 
elliptic flow at intermediate energies may be particulary useful in extracting the information 
on the nuclear effective interaction. While the elliptic flow at the energies higher than the 
transition energy will be useful to extract the medium correction of nucleon-nucleon cross 
sections because two-body collisions play more important role on collective flow at these 



energieP, |2^. The another aim of this work is to investigate the medium correction of 
nucleon-nucleon cross sections through eUiptic flow in heavy ion coUisions at energies from 
Fermi energy to relativistic energies. 

In this letter, we apply the new version of Improved Quantum Molecular Dynamics model 
(ImQMDOS) to study the excitation function of elliptic flow parameters for ^^"^Am +^^'' Au 
at intermediate energies, and through the comparison between measurement and model 
calculations to extract the information on the effective interaction which related to EoS and 
the medium correction of nucleon-nucleon cross sections. The system size dependence of 
transition energies of elliptic flow from ^^Ni Ni to ^^"^Au +1^'^ Au. will also be studied. 

For the convenience of readers, we first give a brief introduction of the ImQMDOS model. 
The main developments of the ImQMD model compared with the usual IQMD model are 
: introducing 1) the isospin independent and dependent surface energy terms in the en- 
ergy density functional, 2)the constraint on the single particle occupation number, and 3) 
the system size dependent wave packet widthj2^. With the ImQMD model, it is able to 



successfully describe the yields of clusters in intermediate energy heavy ion collisions |2j 
In the ImQMDOS model we introduce the full Skyrme potential energy density functional 
except the spin-orbit term in the local interaction part, which allow us to choose various 
Skyrme interactions which describe the ground states of nuclei and saturated nuclear matter 
similarly well but predict rather different properties away from saturated density. 

In the ImQMDOS model, the nuclear local interaction potential energy density functional 
Vioc{p{r)) reads 

, P^^ I 9sur /-^ \2 I 9sur,iso f-^ /■ -,■,2 

2 Po 7 + 1 Po' 2po Po 

+ ( V + 5P^+' + Cp'/')6' + g/-^^, (2) 

Po 

where p, p„, pp are the nucleon, neutron, and proton density, 5 = (p„ — Pp)/(pn, + Pp) is the 
isospin asymmetry. The first two terms in expression (2) are the iso-scalar bulk energy part, 
the third term is the isospin independent surface energy term, the forth term is the surface 
symmetry energy term and the fifth term is the bulk symmetry energy term. The last term, 
called the pr term, is obtained from the pr term of the Skyrme potential energy density 
functional by applying the Thomas- Fermi approximation to the kinetic energy density r and 
thus the explicit momentum dependence is lost. However, the strength of this term gp^. is 
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rather small compared with other iso-scalar terms. The coefficients in expression (2) are 
therefore directly related to the standard Skyrme interaction parameters as 

a 3 P -^7 

^ = 7^(9ti-5t2-4x2t2)po, (3) 
2 o4 

And the A,B and C in the volume symmetry energy term are also given by the Skyrme 
interaction parameters, 

A = -^(xo + 1/2), B = -|(X3 + 1/2), C = -^{^r/'Qsym, (4) 

where Qsym = StiXi — t2(4 + 5x2). The gpr is determined by 

^?p. = ^(3ti + (5 + 4a;2)t2)(f)^/Vf • (5) 

The to,ti,t2,t3 and xo,xi,X2,X3 in expressions (3)-(5) are the parameters of Skyrme force. 
In the calculations performed in this work, for the bulk symmetry potential energy density, 
we only take form, which is corresponding to the form of the linear density dependence 
of the symmetry potential energy as people usually use but the symmetry energy coefficient 
is calculated with the full energy density functional given by (2). The symmetry potential 
energy should not play important role in the quantity studied in this work. Furthermore, 
we introduce an explicit momentum dependent term as the same form as that in 2^ , which 
reads 

Umd = 1.57[/r2(l + 5 x 10-''Ap^)fp/po (6) 

as we find the explicit momentum dependent term is important for elliptic flow. This term 
provides an effective mass m*/m=(l + m/pdUmd/dp), which is about 0.75 at the Fermi 
momentum and about 0.95 at relative momentum around 800 AMevj3]. The calculations 
show that without explicit momentum dependent term, the behavior of the calculated exci- 
tation function of elliptic flow are not consistent with that of experiments no matter which 
interaction is adopted. This finding is in agreement with the conclusion obtained in 
The Coulomb interaction potential energy are also introduced. By using present model we 
are able to directly test effective interactions, specifically, the various Skyrme interactions 
characterized by different 'K', and 'm*/m' of EoS by comparing the predictions of different 
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'0, SkM*0, 



Skyrme interactions with measurement of elliptic flow. In this work, SkP[2£ 
SLy7[3r], and Sllljs^ interactions are chosen. The flrst three are with similar incompress- 
ibility K^o ~ 200-229 MeV but with different m*/m, the last one is with ~ 354 MeV. 
Table 1 gives the parameters in the energy density functional(2) and the properties of sat- 
urated nuclear matter for the Skyrme interactions employed in this work. 

TABLE I: parameters in the ImQMD model and the properties of saturated nuclear matter for 
Skryme interactions employed in this work 





SkP 


SkM* 


SLy7 


SIII 


a{MeV) 


-356.20 


-317.40 


-293.97 


-122.75 


l3{MeV) 


303.03 


248.96 


215.03 


55.19 


7 


7/6 


7/6 


7/6 


2 


gsur{MeVfm'^) 


19.47 


21.82 


22.64 


18.26 


gsur,iso{MeVfm'^) 


-11.35 


-5.47 


-2.25 


-4.94 


9pr{MeV) 


0.00 


5.92 


9.92 


6.42 


m* /m 


1.00 


0.789 


0.687 


0.763 


Poo(/"^"^) 


0.162 


0.160 


0.158 


0.145 


as{MeV) 


30.66 


30.68 


32.62 


28.78 


KooiMeV) 


200 


216 


229 


354 


In the collision term, the phenomenological density dependent 


in-medium nucleon-nucleon 



cross sections are taken, which reads 

<n = (1 - VP/PoWjr, (7) 



where (jfj^'^ denots the free nucleon-nucleon scattering cross sections '33], which are isospin 
dependent. In the treatment of Pauli- blocking in collision part, neutrons and protons are 
treated separately and two criteria are used as in 



and 

PMock = i-{i-mi-fj), (9) 

where is the phase space distribution function for nucleon i. 
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The fragments are constructed by means of the coalescence model widely used in the 
QMD model calculations in which particles with relative momenta smaller than Pq and 
relative distances smaller than Rq are coalesced into one cluster(here Rq = 3.0/m and Pq = 
250MeK/c are adopted). Fig.l shows the charge distribution of fragments for ^^"^ Au +^^'^ Au 
at Eheam = 60, 150,4:00AMeV at central collisions, respectively. The experiment data(solid 



symbols) are taken from 



The calculation results shown are for forward angles as 



the same as the experimental data. In the calculations, the SkP Skyrme interaction is used. 
One sees from the figure that the calculation results for charge distribution of fragments are 
in good agreement with experimental data. Then, we apply the ImQMDOS model to study 
the excitation function of elliptic flow parameters and try to extract the information on the 
effective interactions and the medium corrections of two-body cross sections. 

Fig. 2 shows the excitation function of elliptic flow parameters at mid-rapidity {\y/ymc'^\ ^ 
0.1) for Z < 2 particles for ^^"^ Au +^^''^ Au collisions at b=5fm( the reduced impact parameter 

1 /3 1 /3 

b/bmax equals to 0.38, and bmax = l-l5{Ap + Aqi )). The calculated elliptic flow is given 
in the rotated reference frame as the same as the experimental data. In the figure, solid 
symbols denote experimental dataj^ 0, Q| and open symbols denote calculation results with 
Skyrme interactions SkP, SkM*, SLy7, and SIII, respectively. Concerning the in-medium 
two-body cross sections, the t] in expression (6) is taken to be 0.2. The general behavior 
of the excitation functions of elliptical flow parameters f 2 calculated with different Skyrme 
interactions are similar, i.e. the elliptic flow evolves from a preferential in-plane (rotational 
like) emission (t>2 > 0) to out-of-plane (squeeze out) emission (t>2 < 0) with increase of ener- 
gies. But the detailed behavior of the results from different Skyrme interactions are rather 
different. One can see from the figure that the transition energies at which the elliptic flow 
parameter(f 2) changes sign from positive to negative are divergent for different Skyrme inter- 
actions. The difference is more than 30 MeV per nucleon among the calculation results with 
Skyrme interactions SkP, SkM*, SLy7, and SIII. The transition from preferential in-plane 
emission to out-of-plane emission is because the mean field which contributes to the forma- 
tion of a rotating compound system becomes less important and the collective expansion 
process based on the nucleon-nucleon scattering starts to be predominant. The competition 
between the mean field and the nucleon-nucleon collisions should strongly depend on the 
effective interaction, which leads to the divergence of the transition energies calculated with 
different Skyrme interactions. Clearly, the harder EoS provides stronger pressure which 
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leads to a stronger out-of plane emission and thus to have a smaller transition energy. The 
transition energies calculated with SkP and SkM* are in agreement with experimental data 
while that with SIII and SLy7 are too small compared with experimental data. To see the 
relation between the elliptic flow and the EoS, in the inset in Fig. 2 we show the pressure 
as a function of density calculated from the potential energy density functional (2) for SkP, 
SkM*, SLy7, SIII interactions, respectively. One can see that the transition energy sensi- 
tively depends on the stiffness of EoS, which depends on both K and m* /m. Thus, the best 
flt to the transition energy of the elliptic flow at intermediate energies provides us with the 
information on the stiffness of EoS. It seems to us that one needs multi-observable in order 
to extract K and m*/m explicitly(also seej^, 16|). 

As energy further increases V2 becomes negative, and it reaches maximal negative value 
around 400 AMeV for SkP and SkM* and 250 AMeV for SLy7 and SIII. The calculations 
with SIII and SLy7 provide stronger pressure at the compression zone compared with SkP 
and SkM*, which makes calculated elliptic flow to reach the maximal negative V2 at lower 
energy for SIII and SLy7. To compare the predictions with 4 Skyrme interactions with 
measurement we flnd that the results with SkP and SkM* are in reasonable agreement with 
experimental data. After reaching the maximal negative elliptic flow, the negative V2 value 
decreases again. It irnplies that the spectator moves faster after the f 2 reaches the maximal 
negative valuej^. In 37 1, the nuclear stopping from 90 AMeV to 1.93 AGeV was measured 
and maximal nuclear stopping was observed around 400 AMeV for ^^'^Au+^^''^Au. It seems to 
us that the energy for reaching the maximal negative elliptic flow parameter is in coincident 
with the energy for reaction reaching the maximal nuclear stopping. It is clear that if the 
reaction system reaches the maximal stopping around certain energies the matter formed in 
the reaction should reach minimal transparency and thus particles are preferentially out of 
plane emission mostly. 

Now let us investigate the influence of the medium correction of nucleon-nucleon cross 
sections on elliptic flow. Fig. 3 shows the excitation functions of elliptic flow parameters cal- 
culated with 7] = 0.2, 0.0, —0.4 in expression (6) by which we effectively study the medium 
correction of nucleon-nucleon cross sections at different nuclear environment as well as the 
relative momentum of the scattering pair. The SkP Skyrme interaction is adopted in the 
calculations of Fig. 3. From the flgure, we see that at energies lower than transition energy 
the difference between the calculation results with rj = 0.2 and rj = 0.0 is small and both 
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give reasonable agreement with experimental data and the difference increases when the 
bombarding energy is higher than transition energy. As energy further increases the neg- 
ative elliptic flow calculated with rj = 0.2 is too weak( i.e. too small negative elliptic flow 
parameter). One needs a smaller rj or even a negative 1]. We find a reasonable agreement 
with experimental results can be obtained for the case at incident energy around 400 AMeV 
when 7] is taken to be about —0.4, i.e. at the energy about 400 AMeV, the in-medium two- 
body cross section extracted is larger than the free cross section. In j38|, |39| it was predicted 
that the behavior of the in-medium elastic nucleon-nucleon cross section at super-normal 
densities as a function of the relative momentum of two colliding nuclei is first suppression 
and then enhancement. It is also predicted that the in-medium elastic nucleon-nucleon cross 
section increases with temperature. If we simply consider the relative momentum of col- 
liding nucleon pair to be roughly equal to the relative momentum of projectile and target, 
and suppose the temperature increases obviously from several tens AMeV to several hun- 
dreds AMeV, the information on the in-medium nucleon-nucleon cross sections extracted 
from the elliptic flow is qualitatively consistent with the prediction of This study 

suggests that the rj in the phenomenological expression of in-medium nucleon-nucleon cross 
section (6) should depend on the reaction energy in order to mimic the medium correction 
of nucleon-nucleon cross sections at different environment. To confirm this finding we make 
similar calculations for the excitation function of nuclear stopping in Au+Au at SIS ener- 
gies and compared with measurement . The information about the medium correction of 
two-body cross section extracted from nuclear stopping is in good agreement with that ob- 
tained from excitation function of elliptic flow in this work. The results concerning nuclear 
stopping will be given in another publications. 

We notice that the calculation results are not in fully agreement with measurement at 
whole energy region. It means that a more self-consistent treatment including the in-medium 
cross section and the mean field, especially, a more self-consistent explicitly momentum 
dependent term is needed but it is still difficult up to now. 

We further carry out the study of the system size dependence for the elliptic flow of 
Z <2 particles in ^^Ni Ni, ^^"^Sn +^^2 ^^^Au +^^^ Au. We find that the transition 
energies for three systems are obviously different. We then make a systematic investigation 
of the system size dependence of the transition energies of elliptic flow at intermediate energy 
regime. Fig. 4 shows the transition energies as a function of masses of combined systems. All 
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reactions calculated are of symmetric reactions, the reduced impact parameters are chosen 
to be 0.38. The SkP Skryme interaction and rj — 0.2 in the phenomenological expression(6) 
are adopted. Prom the figure, one sees that the transition energy decreases with the reaction 
system size increasing. One of the important reason is because the pressure produced by 
Coulomb interaction increases with the system size. We fit this curve with the following 
power law, 

Etran^x{Ap + AT)-^. (10) 

The exponent r is about 0.223. Here the exponent is substantially smaller than the exponent 
of the size dependence of balance energies for directed flow. Presumably, it is because more 
complex effects such as the expansion of the compressed zone and the shadowing effect of 
the colder spectator matter play role in changing the sign of elliptic flow compared with the 
directed flow. 

In summary, we have investigated the elliptic flow in heavy ion collisions at energies from 
several tens AMeV to several hundreds AMeV with the ImQMDOS model. By changing the 
Skyrme interactions we study the influence of the EoS on the elliptic flow, especially on the 
transition energy and the energy that the elliptic flow parameter reaches the maximal neg- 
ative value. We flnd that the SkP and SkM* interactions can better describe the excitation 
function of elliptic flow at intermediate energies. The medium correction of nucleon-nucleon 
cross sections is also studied by changing the parameter rj in expression (6). By fltting 
the experimental excitation function of elliptic flow parameters we obtain the behavior of 
in-medium two nucleon cross sections as function of relative momentum of two colliding 
nucleons. Our study suggests that medium correction (the r] value) in the phenomenological 
expression of in-medium cross sections should depend on the relative momentum of colliding 
pair and the medium density and temperature of nuclear medium. The linear density depen- 
dence of the in-medium nucleon-nucleon cross section in (6) probably be better valid when 
the incident energies lower than lOOMeV per nucleon for HIC with heavy nuclear systems. 
The system size dependence of the transition energies is investigated, which flts a power of 
system size with a exponent of 0.223. 
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Fig . 1 The charge distributions of products in the central coUisions of reactions ^^^Am+^^^Am 
at Ebeam=QO, 150, 400AMeV calculated with the ImQMDOS model, respectively. The 
SkP Skyrme interaction is chosen. The experiment data(solid symbols) are taken from 
jsS,!^^. The calculation results shown are for products at forward angles as the same 
as the experimental data. 

Fig. 2 The excitation functions of elliptic flow parameters at mid-rapidity for Z < 2 par- 
ticles from mid-central collisions of ^^'^Au -1-^^'' Au calculated with SkP, SkM*, SLy7, 
SIII Skyrme interactions,respectively. The calculated results are given in the rotated 
reference frame as the same as the experimental data. The experimental data are 
taken from The inset shows the pressure as a function of density calculated with 
SkP, SkM*, SLy7, SIII Skyrme interactions, respectively. 

Fig . 3 The excitation functions of elliptic flow parameters at mid-rapidity for Z < 2 particles 
from mid-central collisions of ^^"^Au +^^'^ Au with rj = 0.2, 0.0, —0.4 in the phenomeno- 
logical expression of in-medium cross sections(6), respectively. The SkP interaction is 



chosen. The experimental data are taken from 



Fig. 4 The transition energies for elliptic flow at intermediate energies as a function of 
combined system mass. 
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